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Given growing interest in optical-frequency magnetic dipole transitions, we use intermediate cou- 
pling calculations to identify strong magnetic dipole emission lines that are well suited for ex- 
perimental study. The energy levels for all trivalent lanthanide ions in the 4/ n configuration are 
calculated using a detailed free ion Hamiltonian, including electrostatic and spin-orbit terms as 
well as two-body, three-body, spin-spin, spin-other-orbit, and electrostatically correlated spin-orbit 
interactions. These free ion energy levels and eigenstates are then used to calculate the oscillator 
strengths for all ground-state magnetic dipole absorption lines and the spontaneous emission rates 
for all magnetic dipole emission lines including transitions between excited states. A large number 
of strong magnetic dipole transitions are predicted throughout the visible and near-infrared spec- 
trum, including many at longer wavelengths that would be ideal for experimental investigation of 
magnetic light-matter interactions with optical metamaterials and plasmonic antennas. 

PACS numbers: 71.20.Eh,31.15.-p,32.50.+d 



I. INTRODUCTION 

The natural optical-frequency magnetic dipole (MD) 
transitions in trivalent lanthanide ions have attracted 
considerable attention in recent years for their ability to 
interact with the magnetic component of light Al- 
though most light-matter interactions are mediated by 
electric fields through electric dipole (ED) transitions, 
the intra-4/ n optical transitions of the lanthanide series 
are well-known to include strong MD contributions* 1 ^ - — 
Spurred by recent advances in optical metamaterials 
and nanophotonics, researchers have proposed a vari- 
ety of ways to leverage natural MD transitions, e.g. 
as the building blocks for homogeneous negative index 
materials^ and as probes for the local magnetic fieldi^r— 
Experimental studies have also demonstrated how the 
competition between ED and MD processes can be used 
to achieve strong enhancement of MD emission^ and 
to broadly tune emission spectral Numerical investiga- 
tions have shown how the enhanced magnetic field in and 
around metal and dielectric nanostructures can promote 
MD transitions r T 11 ' 21 " — illustrating how near-field en- 
hancements can modify optical selection rules to promote 
higher order (ED forbidden) optical processes.— ^~— 

Recent studies have focused primarily on the visi- 
ble 5 Do 7 F\ MD transition in trivalent Europium 
(Eu 3+ ) and the near-infrared 4 /i3/2 — > A I\5/2 MD transi- 
tion in trivalent Erbium (Er 3+ )»ir— £r— The emphasis on 
these transitions is not surprising, because they have a 
long history of scientific and technological importance. 
The 5 L>o -> 7 Fi MD transition in Eu 3+ near 588 nm 
was first characterized in 1941— and subsequently used 
by Drexhage^ Kunz and Lukosai 4 - in their authorita- 
tive studies of modified spontaneous emission. More re- 
cently, spontaneous emission from the Eu 3+ MD tran- 
sition has served as a reference standard in studies of 
local field effects 1 ^ - — and ligand environments^ The 
Er 3+ 4 /i3/2 — > 4 Ii5/2 transition, emitting near 1.5 pm, is 



widely used for fiber amplifiers in optical telecommunica- 
tion. The ED and MD contributions to this mixed tran- 
sition were investigated as early as 1967 by Weberj 18 i 19 
More recently, Er 3+ has been used to demonstrate modi- 
fications in the local density of optical states^ as well as 
stimulated emission along surface plasmon waveguides^ 

From an experimental perspective though, it would be 
helpful to identify additional MD transitions, especially 
in the near infrared range from 700 — 1000 nm. As com- 
pared to the 588 nm visible transition in Eu 3+ , optical 
nanostructures are much easier to fabricate for longer 
wavelengths, and at longer wavelengths, plasmonic res- 
onances also exhibit higher quality- factors due to lower 
Ohmic losses. In contrast to the 1.5 pm line in Er 3+ , 
transitions at wavelengths shorter than 1000 nm can be 
readily observed with high efficiency using standard sili- 
con photodetectors. 

Table 1 in the canonical paper by Carnall et al— has 
served as a definitive list of MD absorption lines for over 
40 years, and since its publication, this table has been the 
basis for identifying possible MD transitions in various 
trivalent lanthanide ions. However, the use of this table 
to identify MD emission lines for experimental study suf- 
fers from two limitations. First and foremost, the table 
restricts itself to transitions involving ground state en- 
ergy levels, and therefore, does not include potential MD 
transition lines that occur between excited states. Sec- 
ond, Ref. limits the free ion Hamiltonian to only the 
electrostatic and spin-orbit interactions. More accurate 
values of the transition wavelengths, oscillator strengths, 
and spontaneous emission rates can be achieved by in- 
cluding higher order terms. 

In this paper, we explicitly calculate MD transitions 
over all possible excited energy levels in the trivalent lan- 
thanide series. We also implement a more complex model 
for the free ion Hamiltonian, including not only the elec- 
trostatic and spin-orbit interactions but also two-body, 
three-body, spin-spin, spin-other-orbit, and electrostat- 
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ically correlated spin-orbit interactions. This model is 
then used to identify all non-zero MD transitions, high- 
lighting those lines that are most promising for experi- 
mental investigation. Using these results, we then ana- 
lyze the effect of various host materials on the branching 
ratio of specific MD transitions. Additionally, calcula- 
tions of electric quadrupole (EQ) transition rates and 
oscillator strengths have been carried out for complete- 
ness and to differentiate MDs from other higher order 
transitions. 



II. METHOD 

Calculations of MD transitions were made by first con- 
structing a Hamiltonian for all 4/ n electron configura- 
tions. The free ion Hamiltonian used is of the formi^ 9 . 

H FI =H + F k fk + CfA so 

fc=0,2,4,6 

+ aL(L + 1) + /3G(G 2 ) + jG(R 7 ) 
+ £ T%+ £ M h m h 

i=2,3,4,6,7,8 /i=0,2,4 

+ E pf pf- C 1 ) 

/=2,4,6 

This Hamiltonian only considers valence electrons. The 
first term, Hq, denotes the central field Hamiltonian that 
shifts the absolute values of the energy levels but not 
their respective spacings. Given that the scope of this 
paper concerns transitions between levels, and their re- 
spective rates, calculations do not include H - For each 
subsequent term, the leading factor represents a radial 
fit parameter that is determined from experiment, while 
the trailing factor is an angular term that can be calcu- 
lated explicitly from first principles. For instance, F k is 
the radial fit parameter for the electrostatic interaction, 
while fk is the calculated angular portion. The spin-orbit 
interaction is designated by £/ and A so . a, /?, and 7 and 
their respective angular portions L(L + 1), G(G2), and 
G(Rr) are the two-body interaction terms. Three-body 
interactions are denoted by T l and U. A combination 
of both the spin-spin and spin-other-orbit interactions 
are encompassed in the M h and irih terms. P* and pf 
denote the electrostatically correlated spin-orbit interac- 
tion. Note that this Hamiltonian does not include terms 
to account for crystal field effects. Although such terms 
are necessary in the calculations of intra-4/ n ED transi- 
tions, they constitute only a small correction for MD and 
EQ transitions, which are directly allowed in intermedi- 
ate coupling. Therefore, the values calculated here are 
representative quantities that can be used to predict and 
analyze MD transitions in any host material. 

After constructing the angular terms using the meth- 
ods outlined in Appendix A, we then used radial fit pa- 
rameters tabulated in Ref. [30| to construct the full Hamil- 
tonian matrix. This matrix was subsequently diagonal- 
ized to yield the free ion energy levels and the \tp[LS]J) 



eigenstates. L,S, and J represent the total orbital, spin, 
and angular momenta, while we use ip to denote all other 
quantum numbers necessary to define each state. Note 
that we place LS in brackets here to illustrate that they 
are no longer good quantum numbers; eigenstates in in- 
termediate coupling are composed of a linear combina- 
tion of different LS terms with the same total angular 
momentum J. Following standard conventions, we label 
each level in Russell-Saunders ( 2S+1 Lj) notation accord- 
ing to their dominant LS term(s). If no single LS term 
has a fractional contribution greater than 50%, then we 
label the level according to the two largest LS terms. 
Using the complete eigenstates, we perform subsequent 
calculations of oscillator strengths and transition rates 
between all levels. Thus, over the full trivalent lanthanide 
series (4/ 1 — 4/ 13 ), we consider a total of 192,177 possible 
transitions, see Table [I] 

TABLE I. Number of terms, levels, and total transitions for 
given / n configuration. 



Configuration 


if 13 ) 


f 

(Z 12 ) 


f 

(/") 


(/ 10 ) 


/ b 
(/ 9 ) 


f b 
(/ 8 ) 


f 


Number of 
Terms (LS) 


1 


7 


17 


47 


73 


119 


119 


Number of 
Levels (LSJ) 


2 


13 


41 


107 


198 


295 


327 


Number of 
Transitions 


1 


78 


820 


5,671 


19,503 


43,365 


53,301 



III. RESULTS AND DISCUSSION 
A. Magnetic Dipole Absorption Lines 

We first calculate the oscillator strengths for all ground 
state MD absorption lines in the trivalent lanthanide se- 
ries. (The formulas used for this calculation are provided 
in Appendix 151) Our results found 468 non-zero MD ab- 
sorption lines, including 84 transitions between 300 nm 
and 10 (im; the vacuum oscillator strengths, P' MD , of 
these transitions are plotted in Fig. [TJ Table |TT] shows a 
list of the most prominent ground state absorption lines, 
restricted to the energy bounds and minimum oscillator 
strengths used in Table 1 of Carnall et al— 

By comparison, we find 13 additional MD transitions 
that are not listed in Ref. H^- While most of these new 
absorption lines are relatively weak, P' MD < 5 x 10 -10 , 
several exhibit significant MD oscillator strengths, in- 
cluding the 2 F 5 / 2 2 Fj/2 (4.14 urn) transition in Ce 3+ , 
& H$/2 4 ^7/2 (235 nm) transition in Sm 3+ , and the 
7 Fq — s- 5 Di (241 nm) transition in Eu 3+ that have vac- 
uum oscillator strengths of 5.24 x 10 -8 , 1.9 x 10~ 9 and 
2.9 x 10 -9 respectively. As well as finding additional ab- 
sorption lines, these calculations provide a more accurate 
prediction of transition wavelengths. For example, the 
4 ^i5/2 — *■ 4 ^i3/2 transition in Er 3+ is here calculated to 
occur at 1528 nm, closer to the observed 1543 nm center 
wavelength 1 ^ than the 1513 nm value reported in Ref. I20I. 
However, it is worth noting that the oscillator strengths 
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TABLE II. Calculated MD vacuum oscillator strengths for trivalent lanthanidesQ 
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a Only transitions with vacuum MD oscillator strength P' MD > 0.015 arc listed. 
b Italic values shown for comparison arc taken from Table 1 of Ref |20i 

c The MD oscillator strength, Pmd, inside a host material with refractive index n r would be: Pmd = P'md nr 



are not significantly changed by the inclusion of higher 
order terms in the free ion Hamiltonian, as evidenced by 
the side-by-side comparison of P' MD values in Table [TTJ 
As further validation, our values also compare favorably 
with the Hartree-Fock code developed by R.D. Cowan 
and maintained at Los Alamos National Laboratory, 40 
which predicts that the 4 /i5/2 — ► 4 -^i3/2 transition in the 
4/ 11 configuration of Er 3+ should occur at 1495.5 nm 
with an oscillator strength of 31.75x 10~ 8 , which is within 
2% of our calculated value of 31.14xl0 -8 . For reference, 
a tabulated version of the all non-zero MD ground state 
absorption lines between 300 and 10 /zm is provided in 
Table SI of the Supplemental Material 4i 



B. Magnetic Dipole Emission Lines 

Beyond ground state absorption lines, there are MD 
transitions that occur solely between two excited states. 
Some of these excited transitions, such as the 5 Do —> 7 ' F\ 
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FIG. 1. (Color online) Plot of the magnetic dipole ground 
state absorption lines and corresponding MD oscillator 
strengths for all trivalent lanthanide ions between 300 and 
10000 nm. 
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FIG. 2. (Color online) Magnetic dipole emission lines and cor- 
responding vacuum emission rates for all trivalent lanthanidc 
ions between 300 and 1700 nm. Strong emission lines with 
vacuum rates greater than 5 s _1 located above the dashed 
line are listed in Table III. 



transition in Eu 3+ and 5 Z?4 — > 7 F 5 transition in Tb 3+ 
have been identified experimentally. 6,12 However, there 
have been no exhaustive studies of MD emission in all 
trivalent lanthanide ions. Here, we use calculations to 
perform such a search. We proceed to tabulate all non- 
zero MD emission lines between 300 and 1700 nm. A 
total of 1927 non-zero MD emission lines were found 
throughout the lanthanide series. In Tables S2-S13 of 
the Supplemental Material we provide a complete list of 
all such transitions, grouping them by originating excited 
level to allow for a more convenient comparison in future 
experimental studies. 41 A more condensed table of strong 
transitions with vacuum emission rates, A' MD , greater 
than 5 s _1 is shown in Table HTT1 

As shown in Figure [H there are many strong MD tran- 
sitions thoughout the ultraviolet, visible, and near in- 
frared spectra. In addition to transitions which have been 
previously identified through ground state calculations or 
experimental characterization, there are many more MD 
emission lines which could be of practical interest. 

In the ultraviolet spectrum, MD transitions in Er 3+ , 
Gd 3+ , and Tb 3+ are particularly strong. The Q P^/ 2 ~^ 



8 S 7/2 (301 nm) and e P 7/2 
tions in Gd 3 



S7/2 (307 nm) transi- 



have vacuum emission rates of 23.64 and 
2 K 15/2 -> 4 / 15/2 
(366 nm) transition in Er 3+ has a vacuum emission rate 



30.24 s 1 , respectively. Similarly, the 



of 18.20 



Note that these transitions to the 4 I 



15/2 



ground state in Er 3+ and the 8 S 7 / 2 ground state in Gd 3+ 
could have been inferred from the absorption lines dis- 
cussed in the previous section. However, the strong UV 
transitions in Tb 3+ occur between excited states, such as 
the 5 D -> 7 Fi (378 nm) and 5 D 1 -> 7 F 2 (381 nm) which 
have vacuum emission rates of 29.24 and 20.20 s , re- 
spectively. These 5 Dj — >• 7 Fj + i Tb 3+ transitions are the 
higher level analogues to the experimentally character- 
ized 5 D4 — >• 7 i 7 5 (530 nm) excited state transition. 



Throughout the visible spectrum, there are strong MD 
transitions in Eu 3+ , Ho 3+ , and Tb 3+ . Similar to the 
UV transitions in Tb 3+ , many of the visible MD transi- 
tions in Eu 3+ and Tb 3+ are higher level analogues to the 
previously known 5 Dj Fj+i transitions. In Eu 3+ , 
the well-known 5 D -4 7 F 1 (584 nm) transition has a 
calculated vacuum emission rate of 14.37 s . In addi- 
tion to this yellow emission line, there are also higher 
energy blue and green MD transitions in Eu 3+ , includ- 
ing the 5 L»i ->■ 7 F 2 (550 nm), 5 D 2 -> 7 F 3 (505 nm), and 
5 Z?3 — > 7 Fn (460 nm) that have vacuum emission rates 
near 10 s _1 each. Likewise, in addition to the green 
5 £>4 — > 7 'F5 (530 nm) line and higher ultraviolet transi- 
tions, Tb 3+ also has several blue- violet MD transitions, 
such as 5 D 2 -> 7 F 3 (409 nm) and 5 D 3 -> 7 'F 4 (420 nm) 
which have vacuum emission rates greater than 15 s . 
Trivalent Holmium (Ho 3+ ) also exhibits several strong 
blue MD transitions. Two prominent Ho 3+ transitions 
are the 3 K$ — ► 5 /g (483 nm) ground state transition and 
the 3 H§ — > 5 Ij (449 nm) excited state transition, which 
have vacuum emission rates of 18.48 and 24.71 s _1 , re- 
spectively. 

Most interestingly from an experimental perspective, 
there are also many strong MD transitions in the near- 
infrared spectrum. At these longer wavelengths, the de- 
sign and fabrication of metamaterials^ 2 - - — resonant opti- 
cal antennas; 4 ^— photonic crystals ) 53 i 54 and plasmonic 
waveguides^— are more established. Although some 
transitions in this regime originate from excited states 
that would require deep UV excitation, there are a num- 
ber of transitions in Dy 3+ , Er 3+ , Tm 3+ , and Yb 3+ that 
can be pumped at visible or near-IR wavelengths and are 
thus strong candidates for experimental use. These in- 
clude several ground state transitions that could be iden- 
tified from the absorption line calculations in the pre- 
vious section, including the 4 /i5/ 2 — > 4 ^i3/2 (1528 nm) 
transition in Er 3+ , the 3 H 5 — > 
tion in Tm 3+ , and the 2 F b / 2 
tion in Yb 3+ . Here, we calculate the MD vacuum emis- 
sion rates of these transitions to be 10.17, 14.55, and 
16.59 s -1 , respectively. Our calculations also reveal sev- 
eral promising excited state MD transitions. These in- 



3 #6 (1219 nm) transi- 
2 F 7/2 (976 nm) transi- 



clude the 4 F 



9/2 



3 ^11/2 (734 nm) transition in Dy 



3+ 



the 4 G 4 3 H 5 (784 nm) transition in Tm 3+ , and the 
( 2 H, 4 G) 11/2 -4 4 / 13/2 (832 nm) transition in Er 3 + that 
have vacuum emission rates of 11.72, 22.64, and 14.86 
s _1 , respectively. 

Of the seven strong near-infrared lines identified above, 
the four transitions between 700 and 1000 nm are the 
most promising candidates for immediate experimental 
study. Unlike longer wavelength transitions such as the 
1.5 |xm transition in Er 3+ , these MD transitions occur in 
a spectral region where they are still readily observed by 
silicon photodctectors. (For example, back-illuminated 
CCD cameras such as the Pixis 1024B from Princeton 
Instruments exhibit greater than 50% quantum efficiency 
up to 900nm.) Nevertheless, these transitions also occur 
at sufficiently long wavelengths that resonant plasmonic 
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TABLE III. Calculated MD vacuum spontaneous emission rates0 
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10.51 




5 P> 3 


7 P4 


460 


9.02 




b D 2 


7 P 3 


505 


11.58 




5 Bi 


7 P 2 


550 


12.29 




3 Pi 


5 P 2 


583 


16.01 




5 P<o 


? Pl 


584 


14.37 




3 Po 


5 Pl 


700 


24.63 




( 3 I, 3 H) 6 


5 G 6 


776 


5.51 


Gd a + 


°P 5 /2 


*S 7 /2 


301 


23.64 




6 P 7 /2 


S 7 / 2 


307 


30.24 




"P 3 


'P 3 


306 


6.65 




5 P 2 


7 P 2 


.308 


10.07 




6 Pi 


7 Pi 


310 


14.40 




5 Pi 


7 Po 


312 


8.81 




5 G 6 


7 Pe 


370 


24.35 




5 Po 


7 Pi 


378 


29.24 




5 G 5 


7 P 5 


381 


14.54 




5 Pi 


7 P 2 


381 


20.20 




5 Pi 


7 Po 


392 


8.21 




B G 4 


7 P4 


393 


9.68 




5 G 3 


? P 3 


399 


5.57 




5 P 2 


? P 3 


409 


17.88 




5 P>3 


7 P4 


427 


15.49 




5 P 2 


7 Pi 


430 


7.11 




5 P4 


7 P 5 


530 


14.32 




( 5 P, 3 P) 2 


5 P 3 


766 


17.49 





SLJ 


S L J 


A (nm) 








D 77- 

-^7/2 


533 


5.13 




6 D 


6 7-1 

°P 5 /2 


555 


8.89 




/ 2 T/ 2 T" \ 

( K, L) 15/2 


^15/2 


635 


9.75 




( 4 P, 6 P) 3/2 


6 P5/ 2 


676 


5.94 




P9/ 2 


Pll/ 2 


734 


11.72 




G9/2 
( 2 P, 2 i)l5/ 2 
( 4 P, 4 P) 3/2 


Gn/ 2 
4 M 15/2 


896 


6.64 




1124 


8.33 




4 6 P 3/2 


1170 


5.19 




G9/2 


Gn/2 


1550 


6.21 


Ho^ 


(°K, i, G) 5 


°/ 6 


361 


12.20 




CG, 3 H) 4 


5 h 


411 


9.86 




3 F 4 


5 F 5 


422 


6.73 




3 H 6 


5 ^7 


449 


24.71 




( x G, 3 H)i 


5 /4 


449 


5.05 




3 K S 


5 / 8 


483 


18.48 




3 P 4 


B F 4 


486 


5.71 




3 K 7 


5 / 7 


486 


8.68 




3 P 2 


5 S 2 


511 


6.61 




3 P 4 


5 P 3 


538 


6.47 




( S G 3 H) 5 


% 


543 


8.35 




( 3 F 3 G) 4 


5 Fi 


618 


7.72 




3 F 4 


( 5 G, 3 G) 5 


653 


16.60 




3 Pi 


5 S 2 


601 


6.19 




3 D 3 


R -F-i 

P4 


672 


12.32 




3 is 


3 P"s 


777 


11.42 




Pi 


5 t-i 

p 2 


800 


5.00 




3 T 


3 


811 


5.60 




f3 TT 3 \ 

(■>H,G) 5 


{ G, H)5 


1078 


6.22 




{3 tt 3/~<\ 

(°H,°G) & 




1126 


12.12 




( F, G) 2 

( 5 n, 5 G) 4 


5 TT' 


1270 


6.40 




( G, ffjs 


1438 


6.48 


T7v. 3 + 

rjl 


^P"l5/2 


J 15/ 2 


366 


18.20 




2 f 
A13/2 


J l 3 / 2 


384 


5.25 




( 2 H, 2 G) 9/2 


4 r 
J ll/ 2 


392 


5.34 




^11/2 


4 7 

J l 3 / 2 


529 


12.36 




2 P>5/2 


4 TP 
r 7/2 


583 


5.05 




2 Pl/ 2 


4 c 
J3/2 

4 17 
F 7 /2 


608 


11.31 




2 n 

"5/2 


686 


20.05 




( 2 G, 4 P) 9/2 
/2 p 2 n \ 

I *l u )3/2 

{2 TT 4 

(, (jJll/ 2 

/2 7J 2 ^-Y\ 

( ti, ^19/2 

12 TT 2 r~t\ 

( n, U)g/2 

l2 n 2 n\ 

( P, D)s/2 
( 2 G, 4 F) g/ 2 

/2 p 2 n\ 

( P, P / ) 3 / 2 


4 J 
I 11/2 


733 


8.48 




4 c 
i>3/2 




5.55 




4 r 

J l 3 / 2 


832 


14.86 




( 2 G, 4 P)g/ 2 


843 


11.21 




4 

^11/2 


978 


11.87 




4 ^5/2 


1081 


8.19 




4 17 

^9/2 


1101 


10.35 




4 17 

^ 3 /2 


1111 


8.56 




Gg/2 

4 r 

- / l 3 / 2 


( H, G)n/ 2 


1276 


12.21 




4 T 

J 15/ 2 


1528 


10.17 




4/ ~< 
"7/2 


<J9/ 2 


1533 


6.43 


Tm J + 


C'p/p) 2 


"P 3 


430 


22.93 




( 3 P, 1 D) 2 


3 P 3 


765 


13.97 




X G 4 


3 P5 


784 


22.64 




( 3 p, 1 p) 2 


3 P 2 


808 


9.29 




( 3 p, 1 p) 2 


( 3 P, 1 P) 2 


983 


12.96 




3 P 3 


3 P4 


1155 


10.88 




X G 4 


3 P4 


1167 


5.60 




3 P 5 


3 P 6 


1219 


14.55 


Yb J + 


"F5/2 


"P 7 / 2 


976 


16.59 



(*G,*P) 5/2 

G7/2 
( 4 G, 4 P) 5 /2 



>H 7 



4 7/2 

Gn/2 
( H, G)g/ 2 

G 7 /2 
( 4 P, 4 G)ii /2 

( 4 p, 4 p) 3/2 

Gg/2 

D7/2 

Gg/2 
4 p 
*l/2 

( 4 G, 4 P) 7/2 

4 Pl/2 

( 4 P, 4 G) 5/2 
4 j 

1 ll/2 
4 J 

J l 3 / 2 
4 T 

!l5/2 

( 4 P, 4 G) 5/2 
( 4 P, 4 G) 5/2 

Gn/2 
^3/2 

( 4 p, 4 p)i /2 



9/2 
5/2 



"5/2 

^11/2 



°P 7 /2 
r 11/2 
6 P 5 /2 
^ 9/2 



'H 



' 9/2 
11/2 
*' 3 / 2 

9/ 2 
P l/ 2 
F7/2 
11/2 
13/2 



15/2 
5/2 



'H 



f 7 /2 
l 3 /2 
P/2 
^3/2 



334 
347 
348 
360 
361 
302 
306 
375 
376 
386 
400 
410 
412 
415 
421 
428 
436 
440 
441 
458 
471 
493 
495 
530 



5.71 
8.28 
5.58 

12.78 

15.99 
6.15 
6.84 
5.74 
6.62 
8.45 
9.44 

11.35 
9.28 
8.77 
6.67 
6.72 
5.71 
9.99 

18.83 
5.09 
8.11 

19.49 
7.27 

10.38 



Only transitions between 300 - 1700 nm with vacuum MD spontaneus emission rate A' MD > 5 s 1 are listed. 
The MD spontaneous emission rate, Amd, inside a host material with refractive index n r would be: Amd = A'^ 



and nanophotonic structures can be readily fabricated. 

For experimental studies, it will also be important to 
select appropriate host materials to maximize MD emis- 
sion. In particular, to enhance the MD contribution to 
mixed transitions, it will be helpful for lanthanide ions 
to be substitutional^ doped into centrosymmetric sites 



where ED transitions are strictly forbidden. Table IIVI 
shows the calculated MD branching ratios for the Yb 3+ 
2 -p5/2 2 Fi 11 (976 nm) transition in different host ma- 
terials. These calculations were performed by comparing 
the total decay rate (r tota ; = 1/r), as inferred from ex- 
perimental lifetime data in the literature^ with the MD 
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TABLE IV. Calculated MD Branching Ratios for Yb 3+ Transition in Different Host Materials 



Host 


Measured Lifetime a 


Refractive Index b 


MD Emission Rate 


MD Branching Ratio 




r (ms) 


n r 


A M d (s -1 ) 


d 


LiYF 4 


2.16 


1.455 


51.10 


11.0% 


LaF 3 


2.22 


1.597 


67.57 


15.0% 


SrF 2 


9.72 


1.438 


49.33 


48.0% 


BaF 2 


8.2 


1.473 


53.02 


43.5% 


KCaF 3 


2.7 


1.378 


43.41 


11.7% 


xv I 3 r i u 


2 08 


1 5 


55 99 


IX. 070 


Rb 2 NaYF 6 


10.84 


1.403 


45.82 


49.7% 


BaY 2 F 8 


2.04 


1.521 


58.38 


11.9% 


Y 2 Si0 5 


1.04 


1.79 


95.15 


9.9% 


Y 3 Al 5 Oi 2 


1.08 


1.82 


100.0 


10.8% 


YAI0 3 


0.72 


1.956 


124.2 


8.9% 


Ca 5 (P0 4 ) 3 F 


1.08 


1.63 


71.85 


7.8% 


LuP0 4 


0.83 


1.83 (est.) 


101.7 


8.4% 


LiY0 2 


1.13 


1.82 (est.) 


100.0 


11.3% 


ScB0 3 


4.8 


1.84 


103.3 


49.6% 



a From Table III in Ref. |62 
b From Table II in Ref. lejl 



spontaneous emission rates {Amd = A' MD n 3 )^~— pre- 
dicted from the vacuum rates in Table EH 63 The MD 
branching ratio is thus defined as: Pmd = Amd total- 
Note that the MD branching ratio for this Yb 3+ transi- 
tion varies significantly in different host materials. In 
centrosymmetric hosts such as SrF2, Rb 2 NaYFg, and 
SCBO3, it is possible to have ~50% of all decay pro- 
cesses result in MD emission. In more common materials, 
such as yttrium aluminum garnet (YAG, Y 3 Al 5 0i2), MD 
emission still accounts for ~10% of all decay processes. 

The relatively simple two-level energy structure of 
Yb 3+ means that MD emission can naturally account 
for a significant contribution to the overall decay. Other, 
more complex energy level structures, such as in Dy 3+ 
and Tm 3+ , mean that there are more decay paths 
from any particular excited state. These transitions 
are thus interesting candidates for enhancing MD emis- 
sion. For instance, the lifetime of the 4 Fg/2 excited 
level in Dy 3+ ranges from 300 ps in LiNbOg 64 to 2.36 
ms in Y 3 Sc 2 Ga30i2 (YSGG)^ leading to respective 
branching ratios of 0.35% and 2.77% for the associated 
4 ^9/2 — > 6 -F\i/2 MD transition. Similar branching ratios 
were found by analyzing the 1 G*4 — > 3 if 5 transition in 
Tm 3+ . 66-68 



C. Electric Quadrupole Calculations 

In the multipolar expansion of light-matter interac- 
tions, MD terms are generally included in the same order 
as EQ terms, because they both scale with spatial deriva- 
tives of the electric field. Thus, a common question is 
to what extent EQ transitions compete with MD tran- 
sitions. For completeness, we have calculated the oscil- 
lators strengths for all EQ ground state absorption lines 
and the spontaneous emission rates for all EQ emission 
lines. The EQ oscillator strengths and transition rates 
were found to be significantly smaller than those for MD 
transitions. 




500 



1 000 2000 
Wavelength, X (nm) 



5000 



FIG. 3. (Color online) Plot of all electric quadruple 
ground state absorption lines and corresponding EQ oscil- 
lator strengths for all trivalent lanthanide ions between 300 
and 10000 nm. 



^ 0.15 

o 

< 

ru 

^ 0.10 

c 

.0 

'(/) 

E 

LU 

E 0.05 



♦ Ce 3 * 


OYb 3 * 


#Pr 3 * 


O Tm 3 * 


A Nd 3 * 


A Er 3 * 


■ Pm 3 * 


□ Ho 3 * 


T Sm 3 * 


v Dy 3 * 


* Eu 3 * 


OTb 3 * 


★ Gd 3 * 



□ 





300 



▲ 



500 700 

Wavelength (nm) 



1000 



1500 



FIG. 4. (Color online) Plot of all electric quadrupole emis- 
sion lines and corresponding EQ vacuum emission rates for 
all trivalent lanthanide ions between 300 and 1700 nm. 
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The strongest EQ transition was the ( 5 -D, 5 P)% — > 5 D Q 
transition in Eu 3+ with a vacuum emission rate of 0.17 
s _1 . While the emission rate for EQ transitions scales 
with n^., this rate is approximately 30 times weaker than 
the weakest MD transition presented in Table Most 
transitions mediated by EQ interactions have an emis- 
sion rate on the order of 0.01 s^ 1 and would thus require 
significant enhancement to even be observed. Figures [3] 
and 0] show the vacuum oscillator strengths and emis- 
sion rates, respectively, for EQ absorption lines and EQ 
emission lines. A complete tabulation of all 236 EQ ab- 
sorption lines (Table S14) and all 3079 EQ emission lines 
(Tables S15-S25) between 300 and 1700 nm is provided 
in the Supplemental Material.— These calculations con- 
firm that EQ transitions in trivalent lanthanide ions are 
negligible in comparison to the MD transitions calculated 
above. 



IV. CONCLUSION 

Using a detailed free ion Hamilitonian, we have cal- 
culated all non-zero MD ground state absorption lines 
and corresponding oscillator strengths throughout the 
full trivalent lanthanide series. These values are well doc- 
umented in the literature, and we observed good agree- 
ment between our results and those found in Ref. Ho- 
using this detailed Hamiltonian, we then calculated all 
non-zero MD and EQ emission lines and their respective 
emission rates for all trivalent lanthanide ions. Although 
the EQ emission rates were found to be negligible, our 
calculations revealed vastly more MD emission lines than 
previously identified by ground state calculations or ex- 
perimental investigation. 

In the specific spectral range from 300 - 1700 nm, 
we identified 1927 MD transitions, including 117 lines 
with vacuum spontaneous emission rates A'md > 5 s~ 4 . 
Of these transitions, four were identified as the most 
promising for experimental exploration: 4 Fg/2 — > 6 -F\i/2 
(734 nm) in Dy 3 +, 1 G 4 -> 3 H 5 (784 nm) in Tm 3 +, 
( 2 tf, 4 G) 11/2 -> 4 /i 3 /2 (832 nm) in Er 3 +, and 2 F 5/2 -> 
2 Ft 12 (976 nm) in Yb 3+ . These near-IR transitions oc- 
cur at wavelengths for which resonant devices are easily 
fabricated, yet still emit within the detection range of 
silicon photodetectors. 

We subsequently demonstrated how free ion calcula- 
tions can be used to analyze and predict MD emission 
within a range of host materials. We compared the calcu- 
lated emission rates with experimental lifetime data from 
the literature to approximate MD branching ratios, and 
for the specific case of the 2 F 5 / 2 excited level in Yb 3+ , 
showed how MD emission can account for up to ~50% 
of all decay processes. These calculations highlighted 
the importance of selecting appropriate hosts, especially 
those with high centrosymmetry and refractive indices, 
to maximize MD contributions. 

These results and the associated tables in the Supple- 
mental Material 4 ^ can thus be used to guide the study of 



magnetic light- matter interactions in trivalent lanthanide 
ions. Beyond the well-known MD emission lines in Eu 3+ 
and Er 3+ , there are many permutations of ions and hosts 
in which MD emission can likely be observed. While fur- 
ther study is needed to find the most practical combi- 
nations, these comprehensive calculations provide a solid 
foundation from which to begin this search, and they 
provide a firm set of numbers with which to analyze fu- 
ture experimental data. The tabulated values may also 
be helpful in studying the potential role of MD transi- 
tions in more complex processes such as upconversion 69 
and quantum cutting J- These same calculations can also 
help focus the design of optical structures to enhance 
MD emission. For example, emission wavelengths, tran- 
sition rates, and branching ratios can be used as the 
starting point for simulating the effects of optical anten- 
nas and metamaterials on MD transitions. Combining 
these quantum-mechanical calculations with experimen- 
tal measurements and electromagnetic simulations can 
expand the toolkit with which to access the naturally 
occurring MD transitions of lanthanide ions. 

ACKNOWLEDGMENTS 

We thank S. Cueff, M. Jiang, S. Karaveli, J. Kurvits 
and D. Li for fruitful discussions. Financial support was 
provided by the Air Force Office of Scientific Research 
(PECASE award FA-9550-10-1-0026) and the National 
Science Foundation (CAREER award EECS-0846466). 



Appendix A: Free Ion Hamiltonian 

Closed form expressions of the interaction terms 
used in these calculations are provided below. These ex- 
pressions are well defined through many different publi- 
cations and are provided here for reference purposes. 

1. Coefficients of Fractional Parentage 

When describing a particular term in the Z n config- 
uration, one must realize that there could be multiple 
ways in which to arrive at that term from the Z n_1 config- 
uration. There is an approach to this problem that was 
developed by Giulio Racah)£Lr— which defines the terms 
of the l n configuration in terms of P" 1 . The terms of U 1 ' 1 
are known as the parents of the corresponding daugh- 
ters l n . These coefficients of fractional parentage (CFP) 
need only be calculated once. For this paper, the CFP 
were not calculated directly but an electronic version of 
the tables produced byNielson and Koster 75 was used 
instead^ All subsequent calculations were made using 
these values. The CFP are denoted by (ip{\tp). Due to 
the fact that a particular state might appear in more than 
one configuration, such as in both the 4/ n and 4/ n+2 con- 
figurations, a method to distinguish when a state appears 
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is necessary. This is accomplished by using the seniority 
number, which can take integer values from 1 to 7, indi- 
cating in which 4/ n configuration a state first appears. 



2. Electrostatic Interaction 

The electrostatic interaction occurs between config- 
urations with two or more electrons. This is a result of 
the Coulomb repulsion between the two electrons. It is 
calculated from two single electron wavefunctions. The 
electrostatic interaction is diagonal in both J and S val- 
ues and the matrix elements are found using the following 
expression^! 



(l n i>'L'S\f k \l n i)LS) 



1 



c (k) 



' ( 



1 



2L + 1 



x (l n i>LS U {k) 



lyi's) (i n ipLs 



U (k) l n i)LS) 



n{Al + 2 - n) 

~ **(2i + i)(4j + iy" 

C( k ' is the irreducible tensor defined by Racahj! 2 - and 
f/W is the irreducible tensor tabulated by Nielson and 
Koster. 75 Since we are concerned with / n configurations, 
we used I = 3 for all calculations. Again, we are using 
the notation in which ip represents all other quantum 
numbers that are not specifically mentioned. 



3. Spin-Orbit Interaction 



5. Three-Body Interaction 

The three-body interaction terms are analogous to 
the two-body but exist for only 4/ 3 to 4/ u . The form 
of this operator is£I 



n — 3 



This operator is built up recursively using the values for 
the 4/ 3 states found in Tables 1 and 2 of Ref. l82l . 



6. Spin-Spin Interaction 

The spin-spin interaction is analogous to the spin- 
orbit but is the interaction between the spins of two elec- 
trons. H ss is calculated recursively, using the reduced 
matrix operator T^ 22 ) . y( 22 ) is defined for the 4/ 2 config- 
uration, these defined values then permit the calculation 
for all 4/ n , n > 2, configurations and using the following 
equation^ 3 - 



(rv|^ 22 )|/v) = ^(-i) 



S'+L+J 



X g (V{|V) 



S' V J 
L S 1 



T (22) 



The spin-orbit interaction is, in essence, a dipole- 
dipole interaction. The spin-orbit interaction is diagonal 
in J but not in L or S. We calculated this interaction 
using the following formula: 



(r^'L'S' \A so \ F4>LS) = (-l) J + L + s ' [g, L g j 
x (f n ip'L'S' V {11) f n i>Ls\ 



Here we are using the conventional notation for the 
Racah 6-j symbols and y( n ) is the irreducible tensor tab- 
ulated by. Nielson and Koster— 



4. Two-Body Interaction 

For configurations with two or more valence elec- 
trons (or holes), 4/ 2 to 4/ 12 , two-body interactions are 
used to help correct for the use of single electron wave- 
functions. The first term in this correction was discov- 
ered by. Trees— The other two terms are calculated using 
the Racah numbers and the Casimir operator The 
eigenvalues of the Casimir operator on the groups i?7 and 
Gi can be found in.Wybourne— 



7. Spin-Other-Orbit and Electrostatically 
Correlated Spin-Orbit Interactions 

The spin-other-orbit interaction is an interaction be- 
tween the spin of one electron and the orbit of another. 
It is only valid for 4/ 2 to 4/ 12 configurations. The elec- 
trostatically correlated spin-orbit interaction is a con- 
figuration interaction between the spin of an electron 
in one configuration with the orbit of an electron re- 
siding in a different configuration. These terms were 
grouped together for calculation by Judd, Crosswhite and 
Crosswhite^ 3 . The following form was usedi^ 

(/> + - az 13 \ /V) = 5jj,(-l f +L+J 



Both r' 11 ) and t (11) are reduced matrix operators. These 
reduced matrix operators in addition to the values a and 
013 are defined for the 4/ 2 configuration in Refs. [H and 
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Appendix B: Magnetic Dipole Transitions 
1. Oscillator Strength 



Appendix C: Electric Quadrupole Transitions 



Oscillator Strength 



All MD ground state absorption lines were calcu- 
lated using the following equation^ 



87r 2 m P /n r \ 1 



MD 



3he 2 t 



(t) 2jTi Smd ' 



where Smd is the magnetic dipole transition line 
strength. This line strength is defined as: 

S MD = -^Y,\W\ L + 9e S \^\ 2 > 



where g e is the gyromagnetic ratio of the electron. A 
list of all non-zero absorption lines can be found in the 
Supplemental Material.— 



2. Transition Rates 

All MD emission lines were calculated using the fol- 
lowing equation^ 



.4 



MD 



1 167T 3 /J,o f n r 

2J+1 3h VT 



Smd, 



and all non-zero transitions can be found in the Supple- 
mental Material^ 



All EQ ground state absorption lines were calculated 
using the following equation)^ 



Ieq 



112 n 3 al 
225 a 



~>EQ 



2J + 1' 



where S eq is the electric quadrupole line strength and is 
defined as: 

Seq = (-1) S+L ' +J+ V(2J+1)(2J' + 1) 



J J' 2 
L' L S 



Calculated values for the expectation value of the radial 
wavefunctions for the lanthanide series, (r 2 ), were taken 
from Table 21.8 in Ref. __ A list of all non-zero absorp- 
tion lines can be found in the Supplemental Material^ 
2. Transition Rates 

All EQ emission lines were calculated using the fol- 
lowing equation^ 



.4 



1 



8tt 5 



EQ 



2 J + 1 5/ieo 



(t) Seq - 



There are a total of 3079 non-zero EQ transitions between 
300 and 1700 nm, all such transitions can be found in the 
Supplemental MaterialJi 
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